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The  electrostatic  potential  Vat(r)  in  the  space  around  an  atom  with  nuclear  charge  Z  and 
electronic  density  pat(r)  is  given  by  eq.  (1),  which  invokes  the  spherical  symmetry  of  the  charge 

distribution.1 


Vat(r)  =  Vat(r)  = 


Z  fpat(r')dr' 
r  J  |r'  -r| 


(1) 


For  neutral  atoms,  Vat(r)  decreases  monotonically  with  radial  distance  from  the  nucleus."  Figure  1 
shows  that  lnVat(r)  varies  nearly  linearly  with  radial  distance  in  the  valence  region,  indicating  that 
the  relationship  is  essentially  exponential,  Vat(r)  =  Aexp  (-Br).  In  Figure  1  are  given  lnVat(r) 
equations  for  nine  atoms,  obtained  by  fitting  the  electrostatic  potentials  computed  from  Hartree- 
Fock  atomic  wave  functions.3  These  equations,  all  of  which  have  correlation  coefficients  better 
than  0.997,  cover  the  ranges  r0  <  r  <  3  bohr,  where  r0  was  taken  to  be  the  radius  of  the  inner 
boundary  of  the  valence  shell,  defined  as  the  radial  distance  to  the  outermost  inflection  point  of  the 
average  local  ionization  energy  of  the  atom."1  (Other  criteria  could  also  be  used  for  defining  this 
boundary,3'7  without  significantly  affecting  our  conclusions.) 

The  electrostatic  potential  of  a  molecule  with  electronic  density  pmol(r)  is  given  by  eq.  (2): 


Vmol(r)  =  £ 
A 


r  pmol(r'  )dr' 

J  |r'  -r 


(2) 


ZA  is  the  charge  on  nucleus  A,  located  at  Ra-  Pathak  and  Gadre  have  shown  that  a  molecular 
electrostatic  potential  has  no  local  maximal  (At  nuclear  sites,  it  is  of  course  undefined.)  Along 
the  axis  between  two  bonded  atoms,  therefore,  Vmol(r)  must  decrease  monotonically  with  distance 
from  each  nucleus  and  reach  an  axial  minimum  at  some  point  rm.9  Since  0Vmol/draxiai)rm  =  0. 
an  infinitesimal  charge  5q  placed  at  rm  will  feel  no  electrostatic  force  from  either  direction  along  the 
intemuclear  axis.  We  have  accordingly  explored  the  possibility  that  rm  may  be  a  meaningful  axial 

boundary  point  between  the  two  bonded  atoms. 

We  have  computed  optimized  geometries,  electronic  densities  and  electrostatic  potentials  for 
the  19  molecules  listed  in  Table  I.  The  calculations  were  carried  out  using  a  non-local  density 
functional  procedure  (Gaussian  94, 10  Becke  exchange11  and  Lee,  Yang  and  Parr  correlation 
functionals,12  6-3 1G**  basis  set).  For  the  bends  indicated,  we  found  the  distances  from  each 
nucleus  to  the  minima  of  both  Vmoi(r)  and  (following  Bader  et  a/13)  pmol(r)  along  the  intemuclear 
axis.  As  a  test  to  determine  whether  either  of  these  points  is  an  effective  boundary  indicator,  we 
compare  these  distances  to  literature  values  of  the  covalent  radii  of  the  respective  atoms. 

Overall,  the  distances  to  the  minima  of  Vmoi(r)  are  in  reasonable  agreement  with  the 
covalent  radii,  especially  in  view  of  the  fact  that  the  latter  do  not  refer  to  any  specific  molecule  but 
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rather  were  obtained  from  experimental  structural  studies  of  a  number  of  compounds.14  The 
greatest  deviadons  occur  for  some  of  the  hydrogens.  On  the  other  hand,  the  distances  to  the 
minima  of  pmol(r)  usually  differ  significantly  more  from  the  covalent  radii  and  show  some  major 
discrepancies,  more  than  0.30  A  in  ten  instances.  The  relative  radii  given  by  the  density  minima 
are  often  not  what  would  be  anticipated;  particularly  striking  examples  are  PH3  and  SiH4,  in 
which  the  hydrogens  are  larger  than  the  central  atoms.  The  electron  density  wells  also  tend  to  be 
wider  and  flatter  than  those  of  the  electrostatic  potential.  Table  H  summarizes  the  data  in  Table  I, 
giving  for  each  atom  the  averages  of  the  covalent  radii  obtained  from  the  potential  minima  and  the 
density  minima.  The  former  are  closer  to  the  literature  values,  often  considerably  so,  for  every 
atom. 

These  results  indicate  that  the  location  of  the  minimum  electrostatic  potential  along  the 
intemuclear  axis  does  define  a  meaningful  boundary  between  two  bonded  atoms.  The  less  realistic 
covalent  radii  obtained  from  the  minimum  electronic  density  are  consistent  with  the  unexpectedly 
large  magnitudes  of  the  atomic  charges  found  by  using  this  point  on  the  intemuclear  axis  as  the 
basis  for  partitioning  the  electronic  charge.13,15  (For  a  relevant  critical  analysis,  see  Perrin.16) 

jf  vm°l(r)  at  its  minimum  on  the  intemuclear  axis,  i.e.  Vmol(rm),  can  be  regarded  as 
reflecting  largely  the  electrostatic  potentials  of  the  two  bonded  atoms  A  and  B,  then  within  this 
approximation,  the  point  rm  corresponds  to  the  intersection  of  VA(r)  and  VB(r),  and  their 
contributions  to  Vm®1(rm)  should  be  equal.  (This  does  not  preclude  other  contributions  to 
Vmol(rm).)  One  can  arrive  at  the  same  conclusion  by  noting  that  several  proposed  definitions  of 
electronegativity  express  it  as  an  electrostatic  potential  or  force  created  by  an  atom  at  a  distance 
equal  to  its  covalent  radius.17'19  It  then  follows  from  the  principle  of  electronegativity 
equalization20'22  that  these  potentials  or  forces  should  be  the  same  for  two  bonded  atoms  at  their 
intemuclear  axial  boundary  point. 

To  test  the  proposed  condition  VA(rm)  ~  VB(rm),  we  have  found  the  points  of  intersection 
of  all  of  the  various  pairwise  combinations  of  the  InV(r)  equations  in  Figure  1,  when  the  respective 
nuclei  are  separated  by  the  sum  of  the  literature  covalent  radii  of  the  atoms.  The  objective  was  to 
determine  how  well  the  radial  distances  to  the  points  of  intersection  reproduce  the  literature 
covalent  radii. 

For  those  combinations  in  which  both  atoms  are  in  the  first  row  of  the  periodic  table  (C  - 
F)  or  both  are  in  the  second  row  (Si  -  Cl),  the  differences  between  the  calculated  and  the  literature 
covalent  radii  are  very  small,  the  average  being  0.013  A.  Excellent  agreement  is  also  obtained  for 
the  second  row/third  row  pairs  (Si-Br,  P-Br,  S-Br  and  Cl-Br);  the  average  discrepancy  is  0.015 
A.  The  deviations  are  larger  for  combinations  of  a  first-row  atom  with  one  in  the  second  row  or 
with  bromine;  the  average  difference  is  0.066  A.  Among  these,  the  covalent  radius  of  the  first-row 
atom  is  invariably  overestimated  and  that  of  the  second-row  atom  or  bromine  is  underestimated. 
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This  may  reflect  the  well-known  discontinuities  in  many  properties  between  atoms  of  the  first  and 
subsequent  rows.23'27  Overall,  these  results  do  provide  significant  support  for  the  idea  that 

bonded  atoms  contribute  approximately  equally  to  Vmol(rm). 

The  main  point  that  we  have  tried  to  make  in  this  paper  is  that  the  electrostatic  potential 
minimum  along  the  intemuclear  axis  provides  a  realistic  boundary  point  between  two  bonded 
atoms,  more  so  than  does  the  electronic  density  minimum.  (In  this  context  it  is  relevant  to  recall 
earlier  work  indicating  that  the  locations  of  the  electrostatic  potential  minima  of  monoatomic 
negative  ions  can  be  identified  with  their  ionic  radii.28)  We  have  also  suggested  that  the  two 
bonded  atoms  contribute  approximately  equally  to  the  minimum  potential  along  the  axis.  It  is 
reasonable  to  suppose  that  these  conclusions  are  related  to  the  observation  that  the  covalent  radius 
of  an  atom  is  quite  close  to  the  radial  distance  at  which  its  electrostatic  potential  equals  its  chemical 
potential.29  The  relationship  is  being  investigated. 
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Figure  Caption: 

Figure  1 .  Natural  logarithm  of  calculated  electrostatic  potential,  in  kcal/mole,  plotted  against  radial 
distance,  in  bohr,  for  the  atoms  indicated.  The  equations  describing  these  plots  in  the 
valence  regions  of  the  atoms  are  also  given.  (See  text  for  definition  of  valence  region.) 
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TABLE  I.  Literature  values  of  covalent  radii  compared  to  results  obtained  in  terms  of  axial 
_ minima  of  (a)  electrostatic  potential  and  (b)  electronic  density. _ 


Molecule 

Bond, 

Bond  Len 

gth,  A 

Covalent  radii,  A 

A-B 

Calc. 

Exp.a 

Literature15 

Min.,  Vmo1 

Min.,  pmo1 

ch4 

C-H 

1.099 

1.094 

C: 

0.772 

0.699 

0.686 

H: 

0.30 

0.400 

0.413 

HCN 

C-H 

1.076 

1.065 

C: 

0.772 

0.683 

0.703 

H: 

0.30 

0.393 

0.373 

nh3 

N-H 

1.028 

1.024 

N: 

0.70 

0.656 

0.754 

H: 

0.30 

0.372 

0.274 

H20 

O-H 

0.976 

0.972 

0: 

0.66 

0.623 

0.774 

H: 

0.30 

0.353 

0.202 

SiH4 

Si-H 

1.494 

1.481 

Si: 

1.17 

1.01 

0.724 

H: 

0.30 

0.487 

0.769 

ph3 

P-H 

1.436 

1.427 

P: 

1.10 

0.971 

0.699 

H: 

0.30 

0.465 

0.737 

h2s 

S-H 

1.359 

1.336 

S: 

1.04 

0.921 

0.837 

H: 

0.30 

0.438 

0.521 

ch3nh2 

C-N 

1.480 

1.471 

C: 

N: 

0.772 

0.70 

0.748 

0.732 

0.620 

0.860 

ch3no2 

C-N 

1.520 

1.489 

C: 

N: 

0.772 

0.70 

0.768 

0.752 

0.617 

0.903 

(H3C)20 

C-0 

1.427 

1.410 

C: 

O: 

0.772 

0.66 

0.729 

0.698 

0.512 

0.914 

c2h5oh 

C-0 

1.441 

1.431 

C: 

O: 

0.772 

0.66 

0.738 

0.703 

0.523 

0.918 

h3c-oci 

C-0 

1.439 

1.418 

C: 

0.772 

0.735 

0.513 

0: 

0.66 

0.703 

0.925 

ch3f 

C-F 

1.398 

1.383 

C: 

0.772 

0.724 

0.482 

F: 

0.64 

0.674 

0.916 

ch3ph2 

C-P 

1.896 

1.863 

C: 

P: 

0.772 

1.10 

0.846 

1.05 

1.14 

0.754 

ch3sh 

c-s 

1.857 

1.819 

C: 

0.772 

0.836 

0.875 

S: 

1.04 

1.02 

0.983 

ch2ci2 

C-Cl 

1.814 

1.772 

C: 

Cl: 

0.772 

0.99 

0.829 

0.985 

0.793 

1.02 

(Continued) 


TABLE  I.  Literature  values  of  covalent  radii  compared  to  results  obtained  in  terms  of  axial 
minima  of  (a)  electrostatic  potential  and  (b)  electronic  density  (continued). 


Bond, 

V  /  -  - ~L - 

Bond  Length,  A 

O 

Covalent  radii,  A 

Molecule 

A-B 

Calc. 

Exp.a 

Literature^ 

Min.,  Vmo1 

Min.,  pmo1 

FNNF 

F-N 

1.428 

1.385 

F: 

0.64 

0.691 

0.778 

N: 

0.70 

0.737 

0.648 

FSSF 

F-S 

1.698 

1.635 

F: 

0.64 

0.722 

0.988 

S: 

1.04 

0.976 

0.710 

SiHsCl 

Si-Cl 

2.098 

2.049 

Si: 

1.17 

1.08 

0.759 

Cl: 

0.99 

1.02 

1.34 

aM.  D.  Harmony,  V.  W.  Laurie,  R.  L.  Kuczkowski,  R.  H.  Schwendeman,  D.  A.  Ramsay,  F.  J.  Lovas,  W.  J 
Lafferty,  and  A.  G.  Maki,  J.  Phys.  Chem.  Ref.  Data  8,  619  (1979). 
bJ.  A.  Dean,  Lange's  Handbook  of  Chemistry,  14th  ed.  (McGraw-Hill,  New  York,  1992). 
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TABLE  II.  Literature  covalent  radii  compared  to  averages  of  values  given  in 

Table  I,  correspondingly  to  axial  minima  of  (a)  electrostatic  potential 
and  (b)  electronic  density. 


Atom 

Number  of 

values 

Covalent  radii,  A 

Literature 

Average, 
min.,  Vmo1 

Average, 
min.,  pmo1 

H 

7 

0.30 

0.415 

0.470 

C 

11 

0.772 

0.758 

0.679 

N 

4 

0.70 

0.719 

0.791 

0 

4 

0.66 

0.682 

0.883 

F 

n 

J 

0.64 

0.696 

0.894 

Si 

2 

1.17 

1.05 

0.742 

P 

2 

1.10 

1.01 

0.727 

S 

3 

1.04 

0.97 

0.843 

Cl 

2 

0.99 

1.00 

1.18 

CO  CL  CO  O  u 
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CO 
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o> 

CM 

N 

CD 

LO 

CM 

CM 

O 

CD 

CO 

CD 

CO 

CO 

a > 

CM 

CD 

CD 

CO 

CM 

CM 

c\i 

CO 

i — 

CM 

c\i 

o 

CO 

i 

i 

CD 

CO 

n 

L 

LO 

CD 

<D 

LO 

T— 

CM 

CM 

o 

CM 
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CD 

q 
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CO 
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O 

CM 

CM 

CM 
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CO 
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CO 
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II 

II 

II 

11 

It 

II 

11 
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c 

C 

c 
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Br:  In  V  =  3.275  S  -2.1520r 


